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ABSTRACT 

With Hubble Space Telescope Wide-Field Planetary Camera 2 observations of the Cygnus Loop su- 
pernova remnant, we examine the interaction of an interstellar cloud with the blast wave on physical 
scales of 10^^ cm. The shock front is distorted, revealing both edge-on and face-on views of filaments 
and diffuse emission, similar to those observed on larger scales at lower resolution. We identify indi- 
vidual shocks in the cloud of density n « 15cm~'^ having velocity Vs « 170kms~^. We also find the 
morphologically unusual diffuse Balmer-dominated emission of faster shocks in a lower-density region. 
The obstacle diffracts these shocks, so they propagate at oblique angles with respect to the primary 
blast wave. The intricate network of diffuse and filamentary Ha emission arises during the early stage 
of interaction between the cloud and blast wave, demonstrating that complex shock propagation and 
emission morphology occur before the onset of instabilities that destroy clouds completely. 

Subject headings: ISM: individual(Cygnus Loop) — shock waves — supernova remnants 



1. THE SOUTHEAST KNOT OF THE CYGNUS LOOP 

Supernova remnants and the interstellar medium act 
upon each other reciprocally. The energy of super- 
nova remnants (SNRs) heats and ionizes the interstellar 
medium (ISM), and their blast waves govern mass ex- 
change between the hot, warm, and cool phases of the ISM. 
In turn, the extant ISM determines the evolution of SNR 
blast waves as they propagate through the environment 
it provides. Thus, in order to understand the large-scale 
structure of the ISM, we must also discern the nature of 
shock evolution in inhomogeneous media. 

The Cygnus Loop supernova remnant provides an ideal 
laboratory in which to examine these issues. It is relatively 
unobscured, with E{B - V) = 0.08 (Parker 1967), and 
nearby (440ti[oo Pc; Blair et al. 1999), so 1" corresponds 
to 7 X 10^^ cm. It is a middle-aged supernova remnant, 
not only in terms of absolute lifetime, r — 8000 yr (Leven- 
son et al. 1998, scaling for the revised distance), but more 
importantly because the interaction of the blast wave with 
the ISM dominates its evolution and its appearance at all 
wavelengths. At optical wavelengths, the outstanding fea- 
tures are due to decelerated shocks that propagate through 
dense interstellar clouds. Shocks that are reflected off the 
cloud surfaces propagate back through previously shocked 
material, further heating and compressing it, enhancing 
X-ray emission (Hester & Cox 1986). 

One example of this characteristic interaction is the 
southeast knot. Fesen, Kwitter, & Downes (1992) drew 
attention to this apparently insignificant feature in the 
optical and X-ray emission, which has a small angular ex- 
tent compared to the diameter of the Cygnus Loop. In 



contrast to the SNR as a whole, the optical appearance 
of the southeast knot suggests that it represents only a 
very minor enhancement in the local ISM, or that the in- 
teraction is very young. Fesen et al. (1992) support the 
former interpretation, proposing that this represents the 
late stage of interaction with a small cloud that has been 
completely engulfed by the blast wave and is in the process 
of being shredded by fluid instabilities (Klein, McKee, & 
Colella 1994). X-ray imaging contradicts this interpreta- 
tion, however, because the knot is located at the apex of a 
large-scale (0.5 degree) indentation in the eastern limb as 
traced by very low surface brightness X-ray emission (Gra- 
ham et al. 1995; Levenson, Graham, & Snowden 1999). 
Thus, the obstacle is certainly large, extending at least 5 
pc along the line of sight, and the interaction is at an early 
stage. 

The optical emission is confined to a 2' x 4' region, but 
the apparent insignificance of the southeast knot may be- 
lie its importance. The highest surface brightest optical 
and X-ray emitting regions of the Cygnus Loop are asso- 
ciated with the well-developed reflected and transmitted 
shocks that form in mature cloud-blast-wave interactions 
(e.g., Hester, Raymond, & Blair 1994). The morphology of 
these bright regions is notoriously complex because multi- 
ple shocks are present along the same line of sight. Only 
in a few cases is the geometry unambiguous, such as the 
western limb (Levenson et al. 1996). In addition, the 
development of fluid dynamic instabilities into the non- 
linear regime during the late phase of interaction adds to 
the challenge of interpreting the more prominent regions. 
These difficulties suggest that focussing on very recent in- 
teractions, although intrinsically fainter, may provide use- 



^ Based on observations with the NASA/ESA Hubble Space Telescope, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-2655. 
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ful insights into the sizes, shapes, and density contrasts 
characterizing the clouds in the vicinity of the Cygnus 
Loop. 

In this work, we use the the Hubble Space Telescope 
(HST) to examine the variation in optical line emission on 
scales of O'.'l = 7 x 10^'* cm, assuming a distance of 440 pc. 
Clarifying the geometry of the emitting regions, we trace 
the motion of the blast wave in this cloud interaction and 
quantify the relevant physical processes. We present the 
data in §2, discuss the morphology and physics in §§3 and 
4, respectively, and summarize our conclusions in §5. 
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Fig. 1. — WFPC2 field of view overlaid on soft X-ray image of 
the surrounding field. The bright optical emission is not coincident 
with the associated X-ray enhancement. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained the data from the HST archive at the Space 
Telescope Science Institute. The observations with the 
Wide-Field Planetary Camera 2 (WFPC2) were performed 
1994 November 25 for Program 5774 (Principal Investiga- 
tor J. Hester). Figure 1 illustrates the location and ori- 
entation of these observations on the soft X-ray image of 
the surrounding field obtained with the RO SAT High Res- 
olution Imager (Graham et al. 1995). A total exposure 
of 2200 s in two frames (to facilitate removal of cosmic 
rays) was taken through each of three narrow filters. The 
F502N filter includes [O III] A5007, F656N includes Ha, 
and F673N includes [S II] AA6717 6731. We employed 
standard HST pipeline processing, then used the IRAF^ 
task crrej to remove cosmic rays in the average images. 
We assembled the individual detector images into a sin- 
gle mosaic using IRAF's wmosaic. Figure 2 contains the 
Ha image, with the detectors PCI, WF2, WF3, and WF4 
identified counterclockwise, beginning at the upper left. 
The 0'.'046 pixels of the PCI correspond to 3 x 10^^ cm, 
and the O'.'l WF pixels are equivalent to 7 x 10^"' cm. We 
apply the flux calibration of Holtzman et al. (1995), us- 
ing the 1998 calibration. The system throughput is 0.104, 
0.111, and 0.052 for Ha, [S II], and [O III], respectively 

^ IRAF is distributed by the National Optical Astronomy Observatories. 
Inc., under cooperative agreement with the National Science Foundatio 



(HST Data Handbook 1998). 

3. MORPHOLOGY 

The supernova blast wave moves from west to east across 
the field of view and has recently encountered a cloud of 
denser-than-avcragc interstellar gas. In the HST field, we 
observe the southern section of the interaction, which ex- 
tends out of the field of view to the north for another 4'. 
We combine the three narrow-band data sets in a false- 
color image (Figure 3), with Hq mapped in red, [S II] 
in green, and [O III] in blue. Magenta thus corresponds 
to strong Ha and [O III], while yellow shows where Ha 
and [S II] appear together. Cyan, which would come from 
strong coincident [O III] and [S II] is almost entirely absent 
from the image. 

In Figure 3, color and morphology are strongly corre- 
lated. The shock morphology is distinctly stratified, with 
several characteristic types of filaments dominating the 
composition of the image. Similar to larger-scale obser- 
vations of the Cygnus Loop at lower resolution, the two 
fundamental characteristic morphologies — sharp filaments 
and diffuse emission — correspond to two distinct viewing 
geometries (Hester 1987). The sharp filaments arise in 
wavy sheets viewed edge-on, through tangencies, and in 
general these shocks propagate in the plane of the sky. 
This is the most favorable geometry for unambiguously 
viewing the stratification of the post-shock flow as it cools 
and recombines. The images therefore reveal in turn the 
Balmer-dominated and [O III]-dominated regions behind 
the shock front. In contrast, shocks viewed face-on pro- 
duce diffuse emission. Because larger swept-up column 
densities are required for this diffuse emission to be de- 
tectable, it is more likely to exhibit the characteristics of 
a complete radiative shock, namely strong [S II] emission 
relative to [O III] and Ha. 

Toward the north and west the morphology becomes 
less filamentary, and secondary colors (magenta and yel- 
low), as opposed to primary colors, become more common, 
with the superposition of multiple cooling stages observed 
along the line of sight. Green and yellow are prevalent 
in the western sections of the southern filaments. This is 
most apparent at the heel of the emission, which delimits 
the western extent of the cloud in WF3, the lower right 
quadrant of the image. The bright yellow complex at this 
location traces strong Ha and [S II], which is character- 
istic of shocks that have swept up a substantial column 
{Nh > 10^*cm~^) and formed complete cooling and re- 
combination zones. 

The slowest shocks show up as [S II] (green) only; oc- 
casionally these regions appear as high surface brightness 
knots (e.g., at the bottom edge of WF2, the lower left 
quadrant). In these instances, the primary blast wave is 
strongly decelerated in the dense cloud medium. More 
common is the extended, faint [S II] emission, which is 
widespread across WF2 and WF3. In general, this rel- 
atively smooth [S II] emission occurs by itself, unassoci- 
ated with Ha or other filaments. This low surface bright- 
ness [S II] emission is bounded on the western side by 
a sharp edge that delineates the current location of the 
shock within the cloud. This edge can be traced from the 

operated by the Association of Universities for Research in Astronomy, 




Fig. 2. — Southeast knot in the light of Hq. The image is scaled linearly from (white) to 4 X lO^-'* crgcm^^ s^^ pixel^-*- (black). Shocks 
viewed face-on produce diffuse emission, while sharp filaments are characteristic of edge-on shocks. The upper left quadrant contains the 
smaller PCI detector, with WF2 at the lower left, WF3 at the lower right, and WF4 at the upper right. 



bright green knots at the western limit of WF2, then run- 
ning north and northwest into WF3, around the bright 
heel and to the northeast up into the bright complex in 
WF4 at the upper right. This edge probably represents 
the original and undisturbed surface of the cloud, since 
the interaction is recent and has not yet had a signifi- 
cant dynamical effect on the cloud. In several locations, 
we identify the initial development of the radiative zone, 
finding [S II] emission downstream of Ha filaments (near 
the right center of WF2 and at the bottom of the heel in 
WF3, for example). The [S II] emission is typically offset 
behind the shock front by Of! 5 (3 x 10^^ cm). Unlike the 
sharply-peaked Ha filaments, the [S Il]-emitting region is 
resolved, with flux extending over an arcsecond-scale re- 
gion of the sky behind any distinct portion of the shock 
front. 

The brightest incomplete [O III] filaments and the [S II]- 
emitting shocks in the heel region are clearly physically 
associated. At the center of the field of view the two 



main Ha/[0 III] filaments arc part of a segment of blast 
wave that is propagating to the southeast. Two tangen- 
cies to this surface form the two most prominent incom- 
plete [O III] shocks and their preceding Balmer filaments. 
The filament turns through 45 degrees to form a funnel-like 
cusp with the heel. Tracing the upper filament towards the 
heel, the [O III] emission first merges with the Ha, produc- 
ing magenta in the false-color image. It then disappears 
as it eventually joins the shocks driven into the western 
extremities of the cloud. Similar morphologies are formed 
by the blast-wave-cloud collision farther north in the XA 
region (Hester & Cox 1986). 

Emission from shocks with well-developed cooling and 
recombination zones comprise the northern section of the 
interaction and fill the WF4 field. The absence of large- 
scale filamentary structures implies that these shocks arc 
more nearly face-on and are lighting up the surface of the 
cloud. The sharp filamentary structures in this region 
have [O III] emission, appearing blue and magenta, and 



4 



Levenson and Graham 




Fig. 3. — False-color image of the southeast knot. Ha emission is mapped in red, [S II] in green, and [O III] in blue. The complete radiative 
cooling zone has developed in bright yellow regions, with strong Ha and [S II]. The bright "heel" at the lower right of the image marks the 
distinct cloud edge. Pure red filaments occur at non-radiative shock fronts. In the examples that have swept up sufficient column density, 
[O III] emission appears close behind, or is dominant, without obvious associated Ho filaments, in regions of incomplete cooling. 



can be connected to the Balmer-doniinated and incom- 
plete shocks farther south. These less-decelerated shocks 
appear in projection against this northern field, so the rel- 
ative east-west position does not correspond to absolute 
advancement within the cloud. 

The prominent, sharp-edged, red structures are Balmer- 
dominated filaments and appear as a part of a contin- 
uous, gently rippled sheet, as in the extreme southern 
section of WF2. These filaments are due to fast shocks 
{vs > lOOkms"^) that excite Ha emission by electron 
collisions in pre-shock gas that is predominantly neutral 
and atomic (Chevalier & Raymond 1978; Chevalier, Ray- 
mond, & Kirshner 1980). The excitation is confined to in 
a narrow zone immediately behind the shock front, and 
the resultant face-on Ha surface brightness is low. Thus, 
in these "non-radiative" shocks, the filaments are seen as 
bright, sharp structures when the shock front is close to 



tangency with the line of sight (cf. Hester 1987). Since 
the gas must be neutral to produce these filaments, they 
also mark regions where the gas is being shocked for the 
first time, delineating the outer boundary of the blast 
wave. Gas-dynamic phenomena in which the gas is mul- 
tiply shocked, including refiected shocks, cannot produce 
Balmer-dominated filaments. 

The Balmer filaments with no associated [O III] (blue) 
emission have swept up Nh < lO^^cm"^. These shocks 
has suffered the smallest deceleration. While the general 
trend in this image is for shocks to propagate from west 
(left) to east (right), filaments at skewed angles reveal 
their direction of propagation when they have swept up 
sufficient column for [O III] {Nh > lO^cm-^) or [S II] 
{Nh ^ 10^^ cm~^) emission to be detectable downstream. 

At the center of the field are two good examples of 
patches of [O III] emission that are correlated with Balmer- 
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line filaments. Their appearance is consistent with a vol- 
ume of [O in] emission bounded on the upstream side by 
the Balmer-line-delimitcd shock transition. The shock at 
this location is propagating to the southeast. Farther east 
in the PCI is an amorphous region of pure [O III]. Any as- 
sociated Ha is very faint, suggesting that the direction of 
propagation of the blast wave here is more nearly face-on. 
The [O III]-dominated incomplete shocks extend to the top 
of the image and are interspersed with and project against 
more complex and diverse emission morphologies. 



4. NON-RADIATIVE FILAMENTS AND INCOMPLETE 

SHOCKS 

In these images, the relationship between the Balmer 
filaments and the downstream [O III] emission of the in- 
complete shock distinguishes the physical conditions that 
are present. We examine the filament near the center of 
the WF3 field, which provides a particularly clear exam- 
ple. At this location, a Balmer-line filament boimds an 
incomplete cooling and recombination zone. The ratio of 
[O III] to Ha surface brightness reaches a plateau of about 
6 at a distance 5 x 10^^ cm behind the current location 
of the shock front, which the center of the Ha emission 
defines. This value exceeds the maximum of about 2 that 
can occur in fully radiative shocks, identifying it unam- 
biguously as an incomplete shock. 
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Fig. 4. — Model [O III] line cmissivity as a function of distance 
behind shock front for various shock velocities and n = lcm~^. 
While in the Vg = lOOkms"^ case, the ijrofilc rises gradually and 
reaches a maximum 10^® cm behind the shock front, in the faster 
examples, the brightest emission is sharply peaked farther behind 
the shock. 

The key parameter that determines the location and 
width of the [O III] emission zone is swept-up column den- 
sity behind the shock front. More exactly, the [O III] sur- 
face brightness profile constrains the ambient density and 
shock velocity. In general, [O III] emission rises gradu- 
ally and closer to slower shock fronts, while the profile is 
sharply peaked and offset farther downstream from fast 
shocks, as Figure 4 illustrates. The effect of increasing 
density is to shift the primary [O III] peak closer to the 
shock front. (In all cases, the initial, smaller peak about 
10^^ cm behind the shock front occurs while oxygen is ex- 
cited to higher ionization states, when [O III] is not the 
dominant coolant.) 
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Fig. 5. — Observed and model [O III] emission profiles. We find 
good agreement between the data (crosses) and model profiles hav- 
ing Vs = 170 (solid line) to 190 (dotted line) kms"'^, assuming line 
of sight depths of 9 x lO^'^ and 1 x 10^* cm, respectively. 

We apply updated versions of the models described by 
Raymond (1979) and Cox & Raymond (1985) for shock ve- 
locities ranging from 100 to 190 kms~^, convolving these 
models with the instrumental profile measured from point 
sources in the field of view. The Balmer filament fixes 
the location of the shock transition, and the free param- 
eters are shock velocity, pre-shock density, and depth of 
the surface along the line of sight. While the models in- 
clude magnetic fields, in the hot, post-shock gas these data 
trace, thermal pressure dominates, and the effect of mag- 
netic pressure is negligible. We find good agreement with 
both Vs = 170 and 190kms~^, with density n = 15 and 
40cm~^ respectively, in these two cases (Figure 5). We 
cannot significantly distinguish between these fast shocks 
because the instrumental resolution dominates the emis- 
sion profile, but we prefer Vs = 170kms~^ because the 
location of peak emission of this model better matches 
the observations. The residuals between the model and 
data at distances around 10^^ cm arc likely the result of 
poor subtraction of the local background where the flux 
is extremely low. Thus, we adopt the model parameters 
Vs = 170 km s"^ and n = 15 cm""^. In calculating the total 
flux, we initially assume that the geometry of the emis- 
sion region is a thin sheet that extends 3 x 10^^ cm in 
the line of sight, which is the observed extent of the fila- 
ment across the image. Comparison of the observed and 
predicted intensities implies that the line of sight depth 
of the filament is in fact factors of 3 to 4 times its extent 
in the plane of the sky. This is a natural result because 
our selection bias favors bright filaments, which tend to 
be those with the greatest extent along the line of sight. 
The inferred velocity is robust against uncertainties in the 
assumed distance to the Cygnus Loop, while the derived 
densities are inversely proportional to distance. 

For comparison, we compute the pressure from the X- 
ray observations of the same region. The surface bright- 
ness measured with the i?05^T High Resolution Imager 
corresponds to an emission measure of 70 cm~^ pc. Here 
we assume a temperature T = 2 x 10^ K and solar abun- 
dances with two caveats; the temperature in this small 
region may be somewhat different than the average tem- 
perature, which is weighted by the brightest X-ray emis- 
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sion; and gas-phase depletion and subsequent grain de- 
struction may alter the abundances (Vancura et al. 1994). 
Using the line of sight depth 10^* cm determined above 
and a 20% contribution of metals to the electron den- 
sity, He, we require n,, = 16 cm"''. Thus, pressure P = 
8.7 X 10~^ dyne cm~^ , which is comparable to the ram pres- 
sure that drives the cloud shock, based on the best-fitting 
parameters above. 

These Balmer filaments and their associated incom- 
plete shocks are distinct from those observed around the 
perimeter of the Cygnus Loop. In those cases, the fila- 
ments define the nearly- circular undisturbed blast wave 
and are continuous over scales of 40' (Levenson et al. 
1998). The corresponding shocks propagate through a 
much less dense medium than the southeast cloud pro- 
vides. Sankrit et al. (2000), for example, find n = 2- 
4cm~^ and Vg ~ 170kms~^ in a filament on the northeast 
periphery, based on ultraviolet line fluxes and intensity ra- 
tios. Assuming the same initial blast wave properties at 
the northeast limb and in the southeast knot, we conclude 
that the latter must be a younger interaction. 

The overpressure that has developed behind the south- 
cast shock drives it strongly into the cloud at a pressure 
exceeding that of the primary blast wave, as measured in 
other cloud encounters of the Cygnus Loop (Raymond et 
al. 1988; Hester, Raymond, & Blair 1994). The maximum 
overpressure of a cloud interaction is a factor of 3 once 
steady flow is established (McKec & Cowie 1975), or up 
to a factor of 6 when the blast wave encounters a plane 
of material (Spitzer 1982), which is more similar to this 
very early stage of interaction with a large cloud. The 
greater overpressure of the southeast knot is indicated not 
only in comparison with the northeastern filament, but 
also (and more significantly) when the average blast wave 
pressure, Pbw ~ 5 x 10~^° dyne cm^^, derived from global 
X-ray data (Ku et al. 1984) is considered. The large over- 
pressure at the southeast knot therefore indicates highly 
non-steady flow. A shock having Vg = 170 km s"^ is un- 
stable unless the transverse magnetic field B > 10/xG, but 
because the timescale for the development of secondary 
shocks is long (on the order of 10* years), they do not 
appear yet in this case (Innes 1992). Thus, the transient 
nature of the current conditions is expected. Only later 
will a slower, large-scale, coherent shock arise in a more 
developed interaction, similar to the western edge of the 
Cygnus Loop, for example. 

Multiple shocks along the line of sight confuse the [O III] 
profile in several other locations where Balmer emission 
bounds a region of bright [O III], so we cannot directly 



compare with the models, but we can characterize them 
qualitatively. We identify the Balmer filament at the cen- 
ter of the mosaicked field as a region of lower density, be- 
cause the [O III] emission is broader and offset farther 
downstream. Immediately behind the filament at the west- 
center of the WF2 field, the fully-radiative signature of 
[S II] implies that this is a higher density region. The emis- 
sion at the extreme southeast of the WF2 field consists 
exclusively of Ha, which suggests that the shock propa- 
gates through lower density here, perhaps in the extended 
envelope of the cloud. Furthermore, this region includes 
not only the sharp filaments that characterize the edge-on 
view, but also more diffuse emission where the shock sur- 
face is viewed at oblique angles. Thus, this region must 
be intrinsically bright. The emissivity depends linearly on 
Vg and n (Raymond 1991), so the shock velocity is higher 
here. 

5. CONCLUSIONS 

The HST images of the interaction between the Cygnus 
Loop blast wave and an interstellar cloud reveal emission 
variations on the smallest measurable scales (3 x 10^* cm). 
The blast wave cannot be identified as a single, uniform 
entity, but is broken into a complex of interacting shock 
fronts as it encounters the obstacle. We identify the char- 
acteristic morphology of sharp filaments, where a shock 
front is viewed edge-on, and diffuse emission, where the 
view is face-on, on the sub-arcsecond scales that WFPC2 
probes. The [O III] profile immediately behind the shock 
front that Balmer-dominated filaments define reveals rel- 
atively fast shocks {vg « 170 km s"^) in the high-density 
(n « 15cm~^) cloud medium. Balmer filaments with- 
out associated [O III] or [S II] emission arise in slightly 
lower-density regions behind faster shocks. Because this 
emission is intrinsically bright, we detect the diffuse com- 
ponent, which is viewed obliquely, as well as the more 
common sharp filaments where the shock front is viewed 
edge-on. Exhibiting extensive networks of non-radiative 
shocks, the southeast knot must represent the early stage 
of interaction between the cloud and blast wave. This ex- 
ample thus illustrates that complex shock propagation and 
emission morphology occur before the onset of instabilities 
that destroy clouds completely. 

We thank Ravi Sankrit for computing the shock model 
emission profiles. Support for this work was provided by 
NASA through grant number AR-08005 from the Space 
Telescope Science Institute, which is operated by AURA, 
Inc., under NASA contract NAS 5-26555. 
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